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Abstract 



We perform a fit to the full data set corresponding to 33.0 kt-yr of data of the 
Super-Kamiokande experiment as well as to all other experiments in order to 
compare the two most likely solutions to the atmospheric neutrino anomaly 
in terms of oscillations in the i/„ — ► v T and i/„ — > v s channels. Using state-of- 
the-art atmospheric neutrino fluxes we have determined the allowed regions 
of oscillation parameters for both channels. We find that the Am 2 values for 
the active-sterile oscillations (both for positive and negative Am 2 ) are higher 
than for the i/„ — > v T case, and that the increased Super-Kamiokande sample 
slightly favours — > v T oscillations over oscillations into a sterile species v s , 
— > u s , and disfavours — > u e . We also give the zenith angle distributions 
predicted for the best fit points in each of the possible oscillation channels. 
Finally we compare our determinations of the atmospheric neutrino oscillation 
parameters with the expected sensitivities of future long-baseline experiments 
K2K, MINOS, ICARUS, OPERA and NOE. 
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I. INTRODUCTION 



Atmospheric showers are initiated when primary cosmic rays hit the Earth's atmosphere. 
Secondary mesons produced in this collision, mostly pions and kaons, decay and give rise to 
electron and muon neutrino and anti-neutrinos fluxes Q. There has been a long-standing 
anomaly between the predicted and observed / v e ratio of the atmospheric neutrino fluxes 
0. Although the absolute individual or v e fluxes are only known to within 30% accuracy, 
different authors agree that the / v e ratio is accurate up to a 5% precision. In this resides 
our confidence on the atmospheric neutrino anomaly (ANA), now strengthened by the high 
statistics sample collected at the Super-Kamiokande experiment [[§]. This experiment has 
marked a turning point in the significance of the ANA. The most likely solution of the ANA 
involves neutrino oscillations 0. In principle we can invoke various neutrino oscillation 
channels, involving the conversion of into either v e or v T (active- active transitions) or the 
oscillation of v ^ into a sterile neutrino v s (active-sterile transitions). Previously we have 
reported an exhaustive analysis of all available experimental atmospheric neutrino data for 

— > v e and for — > v T channels [0] but we have not discussed the sterile neutrino case. 
This is especially well-motivated theoretically, since it constitutes one of the simplest ways to 
reconcile [PHTOH the ANA with other puzzles in the neutrino sector such as the solar neutrino 
problem [11,12] as well as the LSND result [13] and the possible need for a few eV mass 
neutrino as the hot dark matter in the Universe [ 14|, 15| . Although stringent limits on the 
existence of sterile states with large mixing to standard neutrinos have been obtained from 



cosmological Big-Bang Nucleosynthesis considerations fig] , more conservative estimations 
claim that the relevant effective number of light neutrino degrees of freedom may still allow 



or even exceed N v = 4 [17|, thus not precluding scenarios such as given in |§. Moreover, 
nucleosynthesis bounds might be evaded due to the possible suppression of active-sterile 
neutrino conversions in the early universe due to the presence of a lepton asymmetry that 



can be generated by the oscillations themselves [18 



The main aim of the present paper is to compare the — > v T and the —>■ u s transitions 
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using the the new sample corresponding to 535 days of the Super-Kamiokande data. Our 
present analysis uses the latest improved calculations of the atmospheric neutrino fluxes as 



a function of zenith angle JT9|, including the muon polarization effect |20| and taking into 
account a variable neutrino production point pl| . We determine the allowed regions for 
active-sterile oscillation parameters for the two possible cases Am 2 > and Am 2 < 0. We 
first analyze the data using only the Super-Kamiokande results up to 535 days. Then we 
also perform a global analysis including the results of the Soudan2, 1MB, Frejus, Nusex and 
Kamiokande experiments. Instead of using directly the double-ratio our analysis relies on the 
separate use of the electron and /i type event numbers and zenith angle distributions, with 
their error correlations. Following ref. [0] we explicitly check the agreement of our theoretical 
predictions with the experimental Monte Carlo, rendering reliability on our method. 

Our main result is that the Am 2 values for the active-sterile oscillations (both for positive 
and negative Am 2 ) are higher than for the — > v T case with maximal or nearly maximal 
mixing. We find that the increased Super-Kamiokande sample slightly favours — > v T 
oscillations in comparison with — > u s , and strongly disfavours — > v e . Notice that this 
in contrast with our previous results in Ref. |7[] were the smaller Super-Kamiokande sample 
did not have enought statistical weight to overcome the best fit of the — > v e channel to 
the Kamiokande multi-GeV data. We find that, even though the high Super-Kamiokande 
statistics plays a major role in the determination of the final allowed region, there is still some 
appreciable weight carried by the other experiments. In the global fit case the preference for 
the — > v T channel over the — > v s or z/ M — > v e channels is marginal, in clear contrast to the 
situation where only the Super-Kamiokande data is taken into account. Finally we compare 
our results of the atmospheric neutrino oscillation parameters with the accelerator neutrino 



oscillation searches such as CDHSW \Wj[, CHORUS-NOMAD pEfl , as well as searches at 



reactor neutrino experiments p7|-[29|]. We re-confirm our previous conclusions |7[] that the 
recent data obtained at long-baseline reactor experiment CHOOZ [^] rules out at the 90 
% CL the — > v e channel as a solution to the ANA. We also compare the results of our fit 
with the expected sensitivities of future long-baseline experiments such as K2K |3(J, MINOS 



BU, ICARUS IB3I and NOE IB3I. 



II. ATMOSPHERIC NEUTRINO OSCILLATION PROBABILITIES 

Here we determine the expected neutrino event number both in the absence and the 
presence of oscillations. First we compute the expected number of /x-like and e-like events, 
N a , a = fi,e for each experiment 

N lt = N ltlt + N efl , N e = N ee + N^ e) (1) 

where 

d 

N a/ 3 = n t T / a Ka(h,cos9u,E u )P a p—-e(Ep)dE u dEpd(cos9u)dh . (2) 

J drij u (i[cos U u ) d-t^s 

and P a p is the oscillation probability of v a — > vp for given values of E v ,cos9 v and h, i.e., 
P a p = P(y a — > vp\ E v , cosOv, h). In the case of no oscillations, the only non-zero elements 
are the diagonal ones, i.e. P aa = 1 for all a. 

Here n t is the number of targets, T is the experiment's running time, E u is the neutrino 
energy and $ Q is the flux of atmospheric neutrinos of type a — fi, e; Ep is the final charged 
lepton energy and s(Ep) is the detection efficiency for such charged lepton; a is the neutrino- 
nucleon interaction cross section, and 9 V is the angle between the vertical direction and the 
incoming neutrinos (cos^=l corresponds to the down-coming neutrinos). In Eq. (0), h is 
the slant distance from the production point to the sea level for a-type neutrinos with energy 
E v and zenith angle 9 V . Finally, K a is the slant distance distribution which is normalized to 
one 
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The neutrino fluxes, in particular in the sub-GeV range, depend on the solar activity. In 
order to take this fact into account we use in Eq. (0) a linear combination of atmospheric neu- 
trino fluxes <&™ ax and $™ m which correspond to the most active Sun (solar maximum) and 
quiet Sun (solar minimum), respectively, with different weights, depending on the running 
period of each experiment 0. 

For definiteness we assume a two-flavour oscillation scenario, in which the oscillates 
into another flavour either — > u e , — > v s or — > v T . The Schroedinger evolution 



equation of the — vx (where X = e, r or s sterile) system in the matter background for 
neutrinos is given by p2| 



where 



dt 



Vx J \ Hfix H x J \ vx 
H tl = V fl + cos 20^ , tf x = y x - cos 20^, 

A 777^ 

^=-—81112^, (4) 



4^ 



(5) 



v T = v, = ^{-\y n) , (6) 

K = 0, (7) 



Here Gf is the Fermi constant, p is the matter density at the Earth, M is the nucleon mass, 
and Y e (Y n ) is the electron (neutron) fraction. We define Am 2 = m\—m\ in such a way that if 
Am 2 > (Am 2 < 0) the neutrino with largest muon-like component is heavier (lighter) than 
the one with largest X-like component. For anti- neutrinos the signs of potentials Vx should 
be reversed. We have used the approximate analytic expression for the matter density profile 



in the Earth obtained in ref. | 35fl . In order to obtain the oscillation probabilities P Q/ g we have 
made a numerical integration of the evolution equation. The probabilities for neutrinos and 
anti-neutrinos are different because the reversal of sign of matter potential. Notice that for 
the — > v T case there is no matter effect while for the — > u s case we have two possibilities 
depending on the sign of Am 2 . For Am 2 > the matter effects enhance neutrino oscillations 
while depress anti-neutrino oscillations, whereas for the other sign (Am 2 < 0) the opposite 
holds. The same occurs also for — > v e . Although in the latter case one can also have two 
possible signs, we have chosen the most usually assumed case where the muon neutrino is 
heavier than the electron neutrino, as it is theoretically more appealing. Notice also that, 



as seen later, the allowed region for this sign is larger than for the opposite, giving the most 
conservative scenario when comparing with the present limits from CHOOZ. 

III. ATMOSPHERIC NEUTRINO DATA FITS 



Here we describe our fit method to determine the atmospheric oscillation parameters 
for the various possible oscillation channels, including matter effects for both — > v e and 
Vfj, — > v s channels. In the latter case we consider both Am 2 signs. Some results for the 
Up — > v T and — > v s channels have already been presented in ref. p3 |. However our 



approach is more complete and systematic, complementing that of in many ways. For 



example ref. |23j uses the double ratio of experimental-to-expected ratio of muon-like to 



electron-like events R^ujRj^R and the upward-going/down-going muon ratio p^] . It is well 



known that the double-ratio is not well-suited from a statistical point of view due to its non- 
Gaussian character | |38|| . Although the Super-Kamiokande statistics is better, we prefer to 



rely on the separate use of the event numbers paying attention to the correlations between the 
muon predictions and electron predictions. Moreover, following ref. J7| we explicitly verify in 
our present reanalysis the agreement of our predictions with the experimental Monte Carlo 
predictions, leading to a good confidence in the reliability of our results. We also use the 
atmospheric neutrino flux calculation of ref. [Jl9| instead of the Honda et al fluxes ref. |36| . 



Last but not least we perform a global analysis of the data, instead of focussing only on 
Super-Kamiokande data, admitting however the main role played by the latter experiment. 
In this case we take especial care in implementing the experimental detection efficiencies 



1 37|, as in ref. 0. 

The steps required in order to generate the allowed regions of oscillation parameters were 
given in ref. 0. As already mentioned we focus on the simplest interpretation of the ANA 
in terms of the neutrino oscillation hypothesis. For definiteness we assume a two-flavour 
oscillation scenario, in which the oscillates into another flavour either — > u e , —>■ v s 
or z/ M -> v T . 

As already mentioned, when combining the results of the experiments we do not make 
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use of the double ratio, R^/ e / Rff/e > Du t instead we treat the e and /z-like data separately, 
taking into account carefully the correlation of errors. Following ref. we define the \ 2 

as 

x 2 = E(^r a - < eon/ ) • + ^wjri • W a - < e ^), (9) 

r,J 

where I and J stand for any combination of the experimental data set and event-type 
considered, i.e, / = (A, a) and J = (B,{3) where, A, B stands for Frejus, Kamiokande sub- 
GeV, 1MB,... and a, (5 = e, /i. In Eq. @ ^theory ^ ^ e predicted number of events calculated 
from Eq. (H) whereas Nf ata is the number of observed events. In Eq. @ (j\ ata and crf heory 
are the error matrices containing the experimental and theoretical errors respectively. They 
can be written as 

a 2 u = a a (A)p a ^A,B)a p (B), (10) 

where p a p(A,B) stands for the correlation between the a-like events in the A-type experi- 
ment and /3-like events in 5-type experiment, whereas cr a (A) and crp(B) are the errors for 
the number of a and /3-like events in A and B experiments, respectively. The dimension 
of the error matrix varies depending on the combination of experiments included in the 
analysis. 



We compute p a /3(A, B) as in ref. |p8[ . A detailed discussion of the errors and correlations 
used in our analysis can be found in Ref. . In our present analysis, we have conservatively 
ascribed a 30% uncertainty to the absolute neutrino flux, in order to generously account for 
the spread of predictions in different neutrino flux calculations. Next we minimize the \ 2 
function in Eq. (|9|) and determine the allowed region in the sin 2 29 — Am 2 plane, for a given 
confidence level, defined as, 

x 2 = x l m + 4.61 (9.21) for 90 (99)% C.L. (11) 

There are some minor changes with respect to the assumed errors quoted in ref. [|7J, 
mainly for the Super-Kamiokande experiment. These are the following ||: 



• The error in the /i/e ratio, arising from the error in the charged current cross sections. 
This is estimated to be 4.3%. 

• The error in the electron event number arising from the uncertainty in the neutral 
current cross section, which is estimated to be 3.0% (4.1% ) for sub-GeV (multi-GeV) 
events, respectively. 

• The electron versus muon mis-identification error for Super Kamiokande multi-GeV 
events. This is estimated to be 1.5%. 

In Fig. [I] we plot the minimum value attained by the x 2 function for fixed Am 2 as sin 2 29 
is varied freely, as a function of Am 2 . Notice that for large Am 2 ^0.1 eV 2 , the \ 2 is nearly 
constant. This happens because in this limit the contribution of the matter potential in 
Eq (|5|) can be neglected with respect to the Am 2 term, so that the matter effect disappears 
and moreover, the oscillation effect is averaged out. In fact one can see that in this range we 
obtain nearly the same x 2 f° r the z/ M — > v T and — > v s cases. For very small Am 2 ^ 10~ 4 
eV 2 , the situation is opposite, namely the matter term dominates and we obtain a better 
fit for the — > v T channel, as can be seen by comparing the — > v T curve of the Super- 
Kamiokande sub-GeV data (dotted curve in the left panel of Fig. [I]) with the solid {y^ — > u s ) 
and dashed (z/ M — > u e ) curves in the left panel of Fig. [l|). For extremely small Am 2 ^ 10~ 4 
eV 2 , values \ 2 is quite large and approaches a constant, independent of oscillation channel, 
as in the no-oscillation case. Since the average energy of Super-Kamiokande multi-GeV 
data is higher than the sub-GeV one, we find that the limiting Am 2 value below which \ 2 
approaches a constant is higher, as seen in the middle panel. Finally, the right panel in 
Fig. |l] is obtained by combining sub and multi-GeV data. 

A last point worth commenting is that for the — > v T case in the sub-GeV sample there 
are two almost degenerate values of Am 2 for which \ 2 attains a minimum. From Fig. [I] 
one sees that the corresponding oscillation parameter values are Am 2 = 1.1 x 10~ 4 eV 2 and 
2.4 x 10 _3 eV 2 , both with maximal mixing. For the multi-GeV case there is just one minimum 
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at 1.4 x 10 3 eV 2 . Finally in the third panel in Fig. p] we can see that by combining the Super- 
Kamiokande sub-GeV and multi-GeV data we have a unique minimum at 1.4 x 10 _3 eV 2 . 

IV. RESULTS FOR THE OSCILLATION PARAMETERS 

The results of our % 2 fit of the Super-Kamiokande sub-GeV and multi-GeV atmospheric 
neutrino data are given in Fig. |2|. In this figure we give the allowed region of oscillation 
parameters at 90 and 99 % CL. The upper left, upper right, lower left and lower right panels 
show respectively the — > u T , — > v e (Am 2 > 0), — > v s (Am 2 < 0) and the — > z/ s 
(Am 2 > 0). The thick solid (thin solid) curves show the sub-GeV region at 90 (99) % CL 
regions and the dashed (dot-dashed) curves the multi-GeV region at 90 (99) % CL regions. 

One can notice that the matter effects are similar for the upper right and lower right 
panels because matter effects enhance the oscillations for neutrinos in both cases. In contrast, 
in the case of — > v s with Am 2 < the enhancement occurs only for anti-neutrinos 
while in this case the effect of matter suppresses the conversion in v^s. Since the yield of 
atmospheric neutrinos is bigger than that of atmospheric anti-neutrinos, clearly the matter 
effect suppresses the overall conversion probability. Therefore we need in this case a larger 
value of the vacuum mixing angle, as can be seen by comparing the left and right lower 
panels in Fig. ||. 

Notice that in all channels where matter effects play a role (all, except the upper left 
panel) the range of acceptable Am 2 is shifted towards larger values, when compared with 
the Up — > v T case. This follows from looking at the relation between mixing in vacuo and in 



matter [p9[ . In fact, whenever the magnitude of the matter potential is much larger than the 
difference between the two energy eigenvalues in vacuum, i.e., when V M ^> Am 2 /E, there is 
a suppression of the mixing inside the Earth, irrespective of the sign of the matter potential 
and/or Am 2 . As a result, there is a lower cut in the allowed Am 2 value, and it lies higher 
than what is obtained in the data fit for the — ► v T channel. For example, let us consider 
the cases of — > v T and — > v s for Am 2 > 0. One can see comparing the upper left and 
lower right panels that the values of Am 2 for z/„ — > v s channel for the sub-GeV sample at 



90% CL are in the range from 2.0 x 10~ 4 eV 2 to 7.5 x 10~ 3 eV 2 , whereas for the — > v T 
channel they are in the range from 5.1 x 10~ 5 eV 2 to 6.3 x 10 _3 eV 2 . 

It is also interesting to analyse the effect of combining the Super-Kamiokande sub-GeV 
and multi-GeV atmospheric neutrino data, given in separate in Fig. |2|. The results corre- 
sponding to this case are presented in Fig.|3|, which also indicates the best fit points. The 
quality of these fits can be better appreciated from Tables | and [EI]. Comparing the values 
of xLin obtained in the present work with those of ref. @ we see that the allowed region 
is relatively stable with respect to the increased Super-Kamiokande statistics. However, in 
contrast to the case for 325.8 days, now the — > v T channel is as good as the —>■ u e , 
when only the sub-GeV sample is included, with a clear Super-Kamiokande preference for 
the — > v T channel. As before, the combined sub-GeV and multi-GeV data prefers the 
v Xi where X = r or sterile, over the — > u e solution. 

The zenith angle distribution changes in the presence of oscillations and it is different 
for each channel. The main information supportive of the oscillation hypothesis comes from 
the Multi-GeV muon data of Super-Kamiokande, shown in the upper right panel of Fig. 
One can see that the zenith angle distribution of the electron data is roughly consistent with 
the no-oscillation hypothesis. To conclude this section we now turn to the predicted zenith 
angle distributions for the various oscillation channels. As an example we take the case 
of the Super-Kamiokande experiment and compare separately the sub-GeV and Multi-GeV 
data with what is predicted in the case of no-oscillation (thick solid histogram) and in all 
oscillation channels for the corresponding best fit points obtained for the combined sub and 
multi-GeV data analysis performed above (all other histograms). This is shown in Fig. |j. 
In the upper left (right) panel we show the Super-Kamiokande sub-GeV (multi-GeV) for 
the muon data and for zenith distributions of the muons expected for the three channels 

— > v T ,v ll — > v e and — ► v s . Notice that since the best fit point for v s occurs at 

sin(26 l ) = 1, the corresponding distributions are independent of the sign of Am 2 . In the 
lower left (right) panel of Fig. |] we give the same information for the electron data. In the 
case of oscillations between two neutrino species the only relevant channel to compare to in 
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this case is the — > v e channel. 

It is worthwhile to see why the —>■ v e channel is bad for the Super-Kamiokande Multi- 
GeV data by looking at the upper right panel in Fig. |j. Clearly the zenith distribution 
predicted in the no oscillation case is symmetrical in the zenith angle very much in dis- 
agreement with the data. In the presence of — > v e oscillations the asymmetry in the 
distribution is much smaller than in the — > v T or — > v s channels, as seen from the 
figure. 

V. ATMOSPHERIC VERSUS ACCELERATOR AND REACTOR EXPERIMENTS 

We now turn to the comparison of the information obtained from the analysis of the 
atmospheric neutrino data presented above with the results from reactor and accelerator 
experiments as well as the sensitivities of future experiments. For this purpose we present the 
results obtained by combining all the experimental atmospheric neutrino data from various 
experiments 0|| . In Fig. |5] we present the combined information obtained from our analysis 
of all atmospheric neutrino data involving vertex-contained events and compare it with the 
constraints from reactor experiments like Krasnoyarsk J27], Bugey and CHOOZ [S5|, and 



the accelerator experiments such as CDHSW |2§, CHORUS and NOMAD @. We also 
include in the same figure the sensitivities that should be attained at the future long-baseline 
experiments now under discussion. 

The upper- left, upper-right, lower-left and lower- right panels of Fig. |5] show respectively 
the global fit of all atmospheric neutrino data for the — > v T , — > v e , — > v a with 
Am 2 < and — > v s with Am 2 > 0. In each of these panels the star denote the best fit 
point of all combined data. 

The first important point is that from the upper-right panel of Fig. |5| one sees that the 



CHOOZ reactor |29| data already exclude completely the allowed region for the — > v e 
channel when all experiments are combined at 90% CL. The situation is different if only 
the combined sub-GeV and multi-GeV Super-Kamiokande are included. In such case the 
region obtained (upper-right panel of Fig. ||) is not completely excluded by CHOOZ at 90% 
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CL. Present accelerator experiments are not very sensitive to low Am 2 due to their short 
baseline. As a result, for all channels other than —>■ u e the present limits on neutrino 
oscillation parameters from CDHSW |25|, CHORUS and NOMAD |§ are fully consistent 
with the region indicated by the atmospheric neutrino analysis. Future long baseline (LBL) 
experiments have been advocated as a way to independently check the ANA. Using different 
tests such long-baseline experiments now planned at KEK (K2K) ||30|| , Fermilab (MINOS) 
|H and CERN ( ICARUS H, NOE [||, and OPERA @) would test the pattern of 



neutrino oscillations well beyond the reach of present experiments. These tests are the 
following: 

1. t appearance searches 

(NC/CC) 

2. The NC/CC ratio, V t, t A < i^A\ e — ■ This is the most sensitive test, namely the ratio 

' (NC/CC) f ar 

of the ratios between the total neutral current events over charged current events in 
a near detector over a far detector. For example the curves labelled NC/CC ratio at 
the upper-left panel of Fig. || delimit the sensitivity regions of these experiments with 
respect to this test. 

3. The muon disappearance, CC near /CCf ar . This test is based on the comparison be- 
tween the number of charged current interactions in a near detector and those measured 
in the far detector. For example, the curves labelled Disappearance at the lower panel 
of Fig. |5| delimit the sensitivity regions of the relevant experiments with respect to this 
test. 

The second test can potentially discriminate between the active and sterile channels, i.e. 

— > v T and — > v a . However it cannot discriminate between — > v s and the no- 
oscillation hypothesis. In contrast, the last test can probe the oscillation hypothesis itself. 

Notice that the sensitivity curves corresponding to the disappearance test labelled as 
KEK-SK Disappearance at the lower panels of Fig. [| are the same for the — > v T and the 
sterile channel since the average energy of KEK-SK is too low to produce a tau-lepton in 
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the far detector []. In contrast the MINOS experiment has a higher average initial neutrino 
energy and it can see the tau's. Although in this case the exclusion curves corresponding 
to the disappearance test are in principle different for the different oscillation channels, in 
practice, however, the sensitivity plot is dominated by the systematic error. As a result 
discriminating between — > v T and — > v s would be unlikely with the Disappearance test 
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In summary we find that, unfortunately, the regions of oscillation parameters obtained 
from the analysis of the atmospheric neutrino data on vertex-contained events cannot be 
fully tested by the LBL experiments, when the Super-Kamiokande data are included in 
the fit for the — > v T channel as can be seen clearly from the upper-left panel of Fig. [5|. 
One might expect that, due to the upward shift of the Am 2 indicated by the fit for the 
sterile case, it would be possible to completely cover the corresponding region of oscillation 
parameters. This is the case for the MINOS disappearance test. But in general since only 
the disappearance test can discriminate against the no-oscillation hypothesis, and this test 
is intrinsically weaker due to systematics, we find that also for the sterile case most of the 
LBL experiments can not completely probe the region of oscillation parameters indicated by 
the atmospheric neutrino analysis. This is so irrespective of the sign of Am 2 : the lower- left 
panel in Fig. || shows the — > v s channel with Am 2 < while the v s case with 

Am 2 > is shown in the lower-right panel. 

VI. DISCUSSION AND CONCLUSIONS 

In this paper we have compared the relative quality of the active-active and active-sterile 
channels as potential explanations of the ANA using the recent 535 days sample of Super- 
Kamiokande as well all available atmospheric neutrino data in the sub GeV and Multi GeV 
range (vertex-contained events). Using the most recent theoretical atmospheric neutrino 



x The KEK-SK Collaboration is making a proposal of an upgrade aimed at seeing the tau's [4C] 
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fluxes we have determined the allowed regions of oscillation parameters for — > vx for all 
possible channels X = e, r, s. We find that the Am 2 values for the active-sterile oscillations 
(both for positive and negative Am 2 ) are higher than for the — > u T case. Moreover 
the increased Super-Kamiokande sample slightly favours — > u T oscillations in comparison 
to z/ M — > z/ s , and disfavours z/^ — * z/ e more strongly. In the global fit of all experiments 
(vertex-contained events) we find a slight preference for the — > z/ T channel over the other 
channels, including — > z/ e . Insofar as the z/^ — > z/ e channel is concerned, there is no great 
change in the presently allowed region of all combined experiments with respect to the earlier 
situation: it is completely excluded by the CHOOZ experiment at the 90 %CL. What is new 
is that the — > v e fit for 535 days is now worse than with the smaller sample and slightly 
disfavoured when in the global fit of atmospheric data aione and more strongly disfavoured 
if Super-Kamiokande data alone are combined. 

We also give the zenith angle distributions predicted for the best fit points for each of 
the possible oscillation channel. Using the zenith angle distribution expected for multi-GeV 
Super-Kamiokande data in the presence of oscillation we compare the relative goodness 
of the three possible oscillation channels. This allows one to understand clearly why the 

— > v T and — > v s channels are much better than the — > v e channel which is in any 
case ruled out by CHOOZ. The main support for the to u T oscillation hypothesis comes 
from the Super-Kamiokande Multi-GeV muon data. 

Finally we compare our determinations of the atmospheric neutrino oscillation parame- 
ters with the expected sensitivities of future long-baseline experiments such as K2K, MINOS, 
ICARUS, NOE and OPERA. We have found that, unfortunately, the regions of oscillation 
parameters obtained from the analysis of the atmospheric neutrino data on vertex-contained 
events cannot be fully probed by most of the LBL experiments. Despite the Am 2 values in- 
dicated by the atmospheric data for the sterile case are higher than for the — > v T channel, 
the need to rely on the disappearance test makes it difficult to completely test the region 
of atmospheric neutrino oscillation parameters at the presently proposed LBL experiments. 
From this point of view a re-design of such experiments would be desirable, as recently 
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suggested by NOE 
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TABLES 



Experiment 




Via -> Vt 


Vy, -> ^ S 


^ v e 


Super-Kam 


Xmin 


6.1 


7.7 


7.3 


sub-GeV 


Am 2 ( 10" 3 eV 2 ) 


0.11 


1.6 


1.2 




sin 2 29 


1.0 


1.0 


0.97 


Super-Kam 


~X,rnin 


6.3 


7.9 


10.8 


multi-GeV 


Am 2 ( 10" 3 eV 2 ) 


1.4 


3.5 


25.4 




sin 2 29 


0.97 


1.0 


0.74 



TABLE I. Minimum value of \ 2 an d the best fit point for Super-Kamiokande sub-GeV and 
Multi-Gev data and for each oscillation channel. Notice that for i/„ — > v s the minimum x 2 is 
practically independent of the sign of Am 2 as the minimum is located at maximum mixing angle. 



Experiment 




v» -»■ v T 


Vy -»■ ^ S 


^ -»■ Ve 




Super-Kam 




13.7 


16.5 


21.9 




Combined 


Am 2 ( 10~ 3 eV 2 ) 


1.4 


2.2 


1.4 






sin 2 29 


1.0 


1.0 


0.96 




All experiments 




49.5 


50.8 


50.8 




Combined 


Am 2 ( 10~ 3 eV 2 ) 


3.0 


3.3 


3.0 






sin 2 29 


1.0 


1.0 


0.99 





TABLE II. Minimum value of \ 2 from for the combined data of Super-Kamiokande only (upper 
part of table) and for all experiments (lower part of table) and the best fit point for each type of 
oscillation channels. 
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FIGURES 



SKSub-GeV SK Mul-GeV SK Combined 




10 5 10 3 10 1 1010 4 10 2 1 1010 4 10 2 110 



Am 2 (eV 2 ) 



FIG. 1. Xmin f° r fixed Am 2 versus Am 2 for each oscillation channel for Super-Kamiokande 
sub-GeV and multi-GeV data, and for the combined sample. Since the minimum is always obtained 
close to maximum mixing the curves for — ► v s for both signs of Am 2 coincide. 
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FIG. 2. Allowed regions of oscillation parameters for Super-Kamiokande for the different 
oscillation channels as labeled in the figure. In each panel, we show the allowed regions for the 
sub-GeV data at 90 (thick solid line) and 99 % CL (thin solid line) and the multi-GeV data at 90 
(dashed line) and 99 % CL (dot-dashed line). 
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FIG. 3. Allowed oscillation parameters of Super-Kamiokande data combined at 90 (thick solid 
line) and 99 % CL (thin solid line) for the different oscillation channels as labeled in the figure. In 
each panel the best fit point is marked by a star. 
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FIG. 4. Angular distribution for Super-Kamiokande electron-like and muon- like sub-GeV and 
multi-GeV events, together with our prediction in the absence of oscillation (thick solid line) as 
well as the predictions for the best fit points in each oscillation channel: — ► v s (thin solid line), 
Vn — ► ^ e (dashed line) and i/^ — ► z/j- (dotted line). The errors displayed in the experimental points 
is only statistical. 
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FIG. 5. Allowed oscillation parameters for all experiments combined at 90 (thick solid line) and 
99 % CL (thin solid line) for each oscillation channel as labeled in the figure. We also display the 
expected sensitivity of the future long-baseline experiments in each channel: MINOS, KEK-SK, 
NOE, ICARUS and OPERA, as well as the present constraints of accelerator and reactor ex- 
periments: CHOOZ, Bugey and Krasnoyarsk for the — > v e channel, CDHSW and CHORUS 
(NOMAD) for — ► u x , where x = r or sterile. The best fit point is marked with a star. 
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